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Recently we reported on the infrared behavior of the Landau gauge gluon and ghost dressing 
functions in SU(3) Wilson lattice gluodynamics with special emphasis on the Gribov problem. 
Here we add an investigation of the spectral properties of the Faddeev-Popov operator at f) = 5.8 
and 6.2 for lattice sizes 12 4 , 16 4 and 24 4 . The larger the volume the more of its eigenvalues are 
found accumulated close to zero. Using the eigenmodes for the spectral representation it turns out 
that for our smallest lattice 0(200) eigenmodes are sufficient to saturate the ghost propagator at 
lowest momentum. We associate exceptionally large values of the ghost propagator to extraordinary 
contributions of low-lying eigenmodes. 
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I. INTRODUCTION 

The infrared suppression of the gluon propagator on 
the one hand and the enhancement of the ghost prop- 
agator at low momentum on the other are closelyre- 
lated to the Gribov- Zwanziger horizon condition P,|3J3 
as well as to the Kugo-Ojima confinement criterion [4|. 
Zwanziger pj has worked out that the continuum be- 
havior of both propagators in Landau gauge are conse- 
quences of restricting gauge fields to the Gribov region 
f2, where the Faddeev-Popov operator is non-negative. 
In general, for a given gauge field a gauge orbit has more 
than one intersection (Gribov copy) within Q. However, 
in the infinite volume limit expectation values taken over 
arbitrary representatives of are predicted to be equal to 
those over the fundamental modular region A, the set of 
gauge fields being absolute maxima of the Landau gauge 
functional defined below. On a finite lattice, however, 
this equality cannot be expected Therefore, the Gri- 
bov ambiguity has to be explored in detail on finite lat- 
tices before drawing any conclusion about the infrared 
behavior of thepropagators mentioned. In previous in- 
vestigations [a, EH the influence of Gribov copies on the 
Landau gauge ghost and gluon propagators was studied 
in detail both for SU(2) and SU(3). Whereas the gluon 
propagator was not seen to be influenced by the choice 
of gauge copies the ghost propagator turned out to be 
clearly copy-dependent in the limit of small momenta. 

In the present letter we are asking, how the singular be- 
havior of the ghost propagator is related to the spectrum 
of the Landau gauge Faddeev-Popov (F-P) operator. It 
is exactly there that the Gribov ambiguity must become 
visible. We will demonstrate that the low-lying eigen- 
values move towards zero as the volume increases and 
that their values are sensitive to the choice of Gribov 
copies. We will also discuss the localization properties 
of the corresponding eigenmodes. The mode expansion 
of the momentum space ghost propagator will be shown 
to converge the slower the higher the momentum is. Re- 
turning to the problem of exceptionally large values of 



the ghost propagator ( "exceptional configurations" ) , dis- 
cussed in [EH- we can show that they are related to 
strong contributions of the lowest eigenmodes. For the 
sake of completeness we remind of a recent investigation 
of the spectrum of the Coulomb gauge F-P operator Q . 

We introduce relevant definitions and notations in 
Sec. [n] Simulation details are given in Sec. IIIII In 
Sec. IIVI the spectral properties of the F-P operator are 
discussed. Configurations with exceptionally large val- 
ues for the ghost propagator are studied in Section^ In 
Section IVII we will draw our conclusions. 



II. DEFINITIONS 

In order to study the ghost propagator using lattice 
simulations one has to fix the gauge for each thermalizcd 
SU(3) gauge field configuration U = {U x ^}. We apply 
the Landau gauge condition which can be implemented 
as a search for a gauge transformation g = {g x } 

U XilA -> 9 U x>fl = g x U x ^ gl +jl , g x e SU (3) 
that maximizes the Landau gauge functional 

X /J — 1 

while keeping the link variables U x ^ fixed. 

The functional Fjj [g] has many different local maxima. 
When the lattice volume V is enlarged for a fixed in- 
verse coupling constant (3 or, alternatively, (3 is decreased, 
more and more of these maxima become accessible by an 
iterative gauge fixing process starting from various ini- 
tial random gauge transformations. The set of gauge 
copies { 9 U} which correspond to different (local) max- 
ima of -F[/[g], where U is kept fixed, are called Gribov 
copies in analogy to the Gribov ambiguity in the con- 
tinuum 3]. All Gribov copies belong to the gauge orbit 
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created by U and satisfy the differential (lattice) Landau 
gauge transversality condition (d IJ 9 A^)(x) = where 

(d^)(x) ee J2 ( 9 A»(x + A/2) - "A^x - A/2)) . (2) 



Here 9 A /1 (x+(i/2) defines the non-Abelian gauge potential 
on the lattice, i.e. 
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tracclcss 



(3) 



In the following, we will drop the label g for convenience, 
i.e. we assume U to satisfy the Landau gauge condition 
such that g = 1 maximizes the functional in Eq. Q rel- 
ative to the neighborhood of the identity. 

The F-P operator is the Hessian of the gauge functional 
Eq. 0} and can be expressed in terms of the (gauge-fixed) 



link variables U x ^ as 



M ab = N A ab 

xy / j x,fi 



A - fi ab A 
°x,y D x,y, °x- 



_ (-tab c 



(4) 



with AH = KcTr[{T a ,T b }(Ux 4l + Ux^)] 
Bf^ = 2-mtTr[T b T a Ux.^} , 



Cab 



= 2-W.tTr[T a T b Ux-^] 



Here T a [a = 1, . . . , 8) denote the generators of the su(3) 
Lie algebra satisfying Tr[T Q T b ] = 6 ab /2. 

For each maximum of the gauge functional -Ft/ [5] the 
corresponding F-P operator has only positive eigenval- 
ues A,- > 0, (i = 1, . . . , 8V — 8), besides of its eight trivial 
zero modes, i.e. the gauge-fixed field configurations U lie 
within the Gribov region f2. We expect that the spectral 
properties of the F-P operator differ for different Gri- 
bov copies. This should have consequences for the ghost 
propagator. We will exploit the spectral representation 
of the inverse of the F-P operator for a given gauge field 
U in terms of its real (ascendent) eigenvalues A^ and its 
(normalized) eigenvectors <pi{x) in coordinate space 
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N 

2 — 1 



(5) 



The eigenvectors are given at each lattice point x as 
8-component color vectors (pi with components 4>f(x). 
They are normalized such that J2 X \4>i( x )\ 2 — 1- Tak- 
ing their Fourier transformed vectors $i(fc) for momenta 
fc M 6 (— Lfj,/2, +L fJi /2] and averaging over a Monte Carlo 
(MC) generated ensemble of gauge field configurations 
we have computed the ghost propagator from truncated 
mode expansions 



G n (q) = (G(k\n)) MC 



(6) 



where 



1 1 

G(k\n) = - 2-$ 4 (fc).^(-fc) (7) 

i— 1 



denotes the contribution of the eigenvalues and eigen- 
modes on a given gauge field configuration. Here the 
vector and scalar product notation refers to the color in- 
dices. The Fourier momenta are related to the physi- 
cal momenta by 

= (2/a)sin(7rfc M /L Al ) 

with a and denoting the lattice spacing and the linear 
lattice extension, respectively. 

If for the whole ensemble of configurations all 
N = 8(V — 1) non-trivial eigenvalues and eigenvectors 
were known, the ghost propagator would be determined 
completely, i.e. G(q) = GN(q) at all momenta. We will 
check the convergence with respect to the order n by 
comparing the truncated propagator G n (q) with the full 
one, G(q), obtained by inverting M for a set of plane 
wave sources (with k ^ 0) orthogonal to the trivial zero- 
modes. For details we refer to reference 0. 

From Eq. Q it is evident that the low-lying eigenval- 
ues and eigenvectors have a dominating impact on the 
ghost propagator. However, for a finite lattice size it is 
not a priori obvious what fraction of the F-P spectrum 
is responsible for the enhancement of the ghost propa- 
gator at the smallest available momenta. Zwanziger has 
argued that the F-P operator has very small eigenval- 
ues I n particular, it should have a high density p(A) 
of eigenvalues per unit Euclicean volume at the Gribov 
horizon (for A > 0). This causes the ghost propagator to 
diverge stronger than 1 /q 2 at q = [l| . We estimate the 
eigenvalue density p at small A by 



P(A) 



h(X, A + AA) 
iVAA 



(8) 



the average number h of eigenvalues per gauge-fixed con- 
figuration within the interval [A, A + AA] divided by the 
bin size AA. For normalization the denominator N = 8V 
has been chosen, since the F-P matrix is a, N x N sparse 
symmetric matrix with N linearly independent eigen- 
states. Note, the trivial zero modes are described by 
a 8<5(A)-peak at A = 0. 



III. SIMULATION DETAILS 

For the purpose of this study we have analyzed pure 
SU(3) gauge configurations thermalized with the stan- 
dard Wilson action at two values of the inverse coupling 
constant (3 = 5.8 and 6.2. Cycles of one heatbath and 
four micro-canonical over-relaxation steps were used for 
thermalization. As lattice sizes we used I2 4 , 16 4 and 24 4 . 
To each thermalized configuration U a set of N cp random 
gauge transformations was assigned. Each has served as 
a starting point for a gauge fixing procedure. We have 
applied standard over-relaxation with over-relaxation pa- 
rameter tuned to w — 1.63. Keeping U fixed this iterative 
procedure generates a sequence of gauge transformations 
with increasing values of the gauge functional (Eq. Q). 



3 



Thus, the final Landau gauge is iteratively approached 
until the stopping criterion in terms of the transversality 
(see Eq. (01) 

maxTr [d^A^d^AUx)] < 1(T 14 (9) 

was fulfilled. Consequently, each random start leads to 
its own local maximum of the gauge functional. How- 
ever, certain extrema of the functional are found multiple 
times. This happened frequently for the 12 4 lattices, but 
rather seldom on larger ones. Note, we used the worst 
local violation of transversality as stopping criteria which 
at a first glance seems to be very conservative. However, 
we have found that the precision of transversality is cru- 
cial for the final precision of the ghost propagator at low 
momentum. 

In order to study the dependence on Gribov copies, for 
each U that copy with largest functional value — among 
all iV cp gauge-fixed ones — was stored, labeled as best 
copy (be). The larger N cp , the bigger the likeliness that 
this copy represents the absolute maximum of the func- 
tional in Eq. QJ, a t least with respect to the observable 
in question. Indeed, with an increasing number N cp of 
inspected gauge copies for each U, the expectation values 
of the ghost and gluon propagator evaluated on be copies 
is found to converge more or less rapidly 5j. The first 
gauge copy was also stored, labeled as first copy (fc). 
Obviously, this copy is as good as any other arbitrarily 
selected one. 

On those ensembles of fc and be gauge-fixed config- 
urations the low-lying eigenvalues A of the F-P opera- 
tor and the corresponding eigenmodes have been sepa- 
rately extracted, where we used the parallelized version 
of the ARPACK package 0, PARPACK. To be specific, the 
200 lowest (non-trival) eigenvalues and their correspond- 
ing eigenfunctions have been calculated at j3 = 6.2 using 
the lattice sizes 12 4 and 16 4 (see Table©. On the 24 4 lat- 
tice, due to restricted amount of computing time, only 50 
eigenvalues and eigenmodes have been extracted at the 
same f3. In addition, 90 eigenvalues have been calculated 
on a 24 4 lattice at [3 — 5.8 providing us with an even 
larger physical volume. This allows us to check whether 
low-lying eigenvalues are shifted towards A — » as the 
physical volume is increased. 



No. 


13 


lattice 


# conf 


# copies 


# eigenvalues 


1 


6.2 


12 4 


150 


20 


200 


2 


6.2 


16 4 


100 


30 


200 


3 


6.2 


24 4 


35 


30 


50 


4 


5.8 


24 4 


25 


40 


90 



TABLE I: Statistics of the data used in our analysis. The last 
column lists the number of eigenvalues extracted separately 
on fc and be copies of U. At f3 — 6.2 the corresponding 
eigenmodes were calculated, too. 
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FIG. 1: The frequency h(X) per configuration of the lowest 
(left panels) and second lowest (right panels) eigenvalue A of 
the Faddeev-Popov operator is shown. Filled boxes represent 
the distribution obtained on be gauge copies, while open ones 
represent those on fc copies. 



IV. THE SPECTRAL PROPERTIES OF THE 
F-P OPERATOR 

A. The spectrum of the low-lying eigenvalues 

Let us first discuss the distributions of the lowest Ai 
and second lowest A2 eigenvalue of the F-P operator. 
These are shown for different volumes in Fig. ^ There 
h(X, A + A) refers to the average number (per config- 
uration) of eigenvalues found in the intervall [A, A + A]. 
Open (full) bars refer to the distribution on f c (be) gauge 
copies. 

From Fig. ^ it is quite obvious that both eigenvalues 
Ai and A2 are shifted to lower values as the physical vol- 
ume is increased. In conjunction the spread of A values is 
decreased. This would be even more obvious, if we had 
shown both distributions as functions of A in physical 
units. It is also visible that the two low-lying eigenval- 
ues A* c (i = 1,2) on fc gauge copies tend to be lower 
than those on be copies. However, this holds only on av- 
erage as can be seen from Fig. [21 There the differences 
X\ c — X\ c of the lowest eigenvalues on f c and be gauge 
copies are shown for different lattice sizes at (3 = 6.2 and 
5.8. It is quite evident that there are few cases where 
A^° < A* c , even though F bc > F tc always holds for the 
gauge functional. 

In addition we have checked how the average values 
(A) of the eigenvalue distributions tends towards zero 
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FIG. 2: The differences A? - \{ c of the lowest F-P eigenvalues 
calculated on be and f c representatives for each gauge con- 
figuration are shown. From left to right the lattice sizes are 
12 4 , 16 4 and 24 4 at (3 = 6.2 and 24 4 at (3 = 5.8. 



as the linear extension aL of the physical volume is in- 
creased. As in our previous study Q we followed Ref. 01 
to fix the lattice spacing a. For j3 = 5.8 and 6.2 we used 
a _1 =1.446 GeV and 2.914 GeV, respectively, using the 
Sommer scale ro = 0.5 fm. 

If the low-lying eigenvalues are supplemented with 
physical units it turns out that the average values of their 
distributions tend towards zero stronger than l/(aL) 2 . In 
fact, using the ansatz 



f(aL) = 



C 



(aL) 



2+s 



(10) 



to fit the data of (Xi)/a 2 for different (aL), a positive e 
is found. The parameter of these fits are given in Ta- 
ble [n] and in Fig. [3] we show the data and the corre- 
sponding fitting functions. There one clearly sees, the 
low- lying eigenvalues not only approach zero, but also 
become closer to each other with increasing aL. We 
have addressed the latter issue by fitting the differences 



a 2 f(aL) 


C 


e 


X 2 /ndf 


<Ai> 


0.120(3) 


0.16(4) 


0.7 


<A 2 > 


0.165(4) 


0.24(5) 


1.8 


(As) 


0.290(1) 


0.45(4) 


3.5 


<A 2 ) - (Ai) 


0.045(2) 


0.47(9) 


0.4 


<A 3 ) - <A 2 > 


0.051(1) 


0.88(8) 


0.2 


<A 4 ) - <A 3 > 


0.033(1) 


0.62(33) 


2.0 


(As) - <A 4 > 


0.037(1) 


0.89(1) 


0.003 



TABLE II: The parameter C and e from fitting either the 
averages (Xi)/a 2 or the differences of adjacent average values 
(Ai+i)/a 2 — (\i)/a 2 of the corresponding eigenvalue distribu- 
tions to the ansatz f(aL) — d/ (aL) 2+ei . 



((Ai+i) — (Xi))/a 2 of adjacent average values using the 
same ansatz Eq. (I10f) . In Table HT1 we give the parameter 
of those fits. 

As mentioned in Sec. |D] the eigenvalue density p(A) 
is of particular interest. We have estimated this quan- 
tity according to Eq. JHJ where the bin sizes have been 
reasonably adjusted for the different volumes. In Fig. 0] 
the estimates are shown for two values of j3. There one 
clearly sees the eigenvalue density close to A = be- 
comes a steeper function of A as the physical volume 
becomes larger. It is remarkable that the increase go- 
ing from (3 = 6.2 to (3 = 5.8 on a 24 4 lattice is even 
larger than going from 12 4 to 24 4 at (3 = 6.2, although in 
both cases the physical volume is increased by a factor 
of about 16. 




l/(oL) [mi" 1 ] 

FIG. 3: The average values (Ai)/a 2 (scaled to physical units) 
of the eigenvalues Ai (i = 1, 2, 5, 20) are shown vs. the inverse 
of the linear lattice extension aL. Only eigenvalues on be 
copies are shown. The lines represent fits to the data using 
the ansatz a -2 (A 4 ) = d/(aL) 2+H . 
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FIG. 4: The eigenvalue density p for be copies as a function 
of A estimated on 12 4 and 16 4 lattices at P = 6.2 and on a 
24 4 lattice for P = 6.2 and 5.8. Bin sizes have been chosen as 
small as possible for each lattice size. The points marks the 
middle of each bin and the lines are to guide the eye. 
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FIG. 5: The relative distribution h of IPR values of the 10 (left), the 10 to 20 (middle) and the 20 to 50 (right) lowest eigenmodes 
are shown. Note there is a logarithmic scale for /i(IPR). Each row corresponds to one pair of f3 and lattice size. Filled boxes 
refer to distributions on be gauge copies, while open ones correspond to f c copies. 



B. Localization properties 

Together with the low-lying eigenvalues the corre- 
sponding eigenvectors (j>(x) have been extracted as well. 
It is interesting to evaluate for each vector the inverse 
participation ratio (IPR). The IPR is defined as 



IPR = V 



E 



4>( X )\ 4 



with V = L 



and is a measure for the localization of an eigenvector. 
It enables us to distinguish between eigenmodes with 
approximately uniformly distributed values of |0(x)| 2 
(IPR w 1 ... 2) and specific ones with a small num- 
ber of sites x having large modulus squared |0(a;)| 2 
(IPR ~ O(100)). Note, the eight (trivial) zero modes 
(A = 0) of the F-P operator give all IPR = 1. 

From Fig. [3] we learn that the majority of eigenvectors 
of the Faddeev-Popov operator is not localized. However, 
some large IPR values have been found associated with 
modes among the 10 lowest non-zero eigenmodes. This 
becomes more likely as the volume is increased. So far we 
have no physical interpretation what causes the stronger 
localization in these rare cases. 



C. What fraction of the F-P spectrum is 
dominating the ghost operator? 

As we have mentioned in Sec. [H] there is an obvious 
way to construct the ghost operator, if all eigenvalues \ 



of the F-P operator and the corresponding eigenvectors 
3>i(fc) in momentum space would be available. However, 
their determination for each configuration is numerically 
too demanding. 

Restricting the sum in Eq. Q to the n lowest eigen- 
values and eigenvectors [n <C N = 8V — 8) , we can figure 
out to what extent these modes, i.e. the corresponding 
estimator G n {q) in Eq. ©, saturate the full ghost propa- 
gator G(q) obtained independently for a set of momenta 
by inverting the F-P matrix on plane waves. See our 
recent study [j| for the data of G(q). 

This saturation is shown in Fig. for the lowest qi 
and the second lowest momentum qi available on differ- 
ent lattice sizes for (3 = 6.2. There the values of G n (q) 
have been divided by the values for the full propagator 
G{q) in order to compare the saturation for different vol- 
umes. Since <&i(fc) has been obtained by a fast Fourier 
transformation of the eigenvector <pi{x), all lattice mo- 
menta k are available. Thus G n (q) refers to the average 
over all k giving raise to the same momentum q. The full 
propagator values G(q) at qi(k) and q2(k), however, refer 
to the averages over lattice momenta k = ([1, 0], 0, 0) and 
to k — (1, 1, 0, 0), respectively. 

Let us consider first the lowest momentum q\. We 
observe from Fig. that the approach to convergence 
differs, albeit slightly, for the three different lattice sizes. 
The relative deficit for n < 50 rises with the lattice vol- 
ume. For n > 100 the rate on a 16 4 lattice is even a bit 
larger than that on a 12 4 lattice. Unfortunately, there 
are no data available for n > 50 on the 24 4 lattice. How- 
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FIG. 6: The ratio of the truncated ghost propagator G n {q) 
(in terms of the n lowest F-P eigenmodes and eigenvalues) to 
the full estimate G(q) (taken from shown as a function of 
n for the lowest (51) and second lowest (52) momentum. The 
inverse coupling is (5 — 6.2 and the lattice size ranges from 
12 4 to 24 4 . The data refer to be copies. 



FIG. 7: Scatter plots of MC time history values of the ghost 
propagator G(k) at k = ([1, 0], 0, 0) vs. the lowest F-P eigen- 
value Ai are shown. The upper panels show data at j3 = 6.2 
on a 12 4 (left) and 16 4 (right) lattice, the lower ones on a 24 4 
lattice at /3 = 6.2 (left) and /3 = 5.8 (right). Crosses refer to 
f c gauge copies and filled circles to be copies. 



ever, for the 12 4 and 16 4 lattices the rates of convergence 
are about the same. For example, taking only 20 eigen- 
modes one is definitely far from saturation (by about 
50%) whereas 150 . . . 200 eigenmodes are sufficient to re- 
produce the ghost propagator within a few percent. In 
other words, the ghost propagator at lowest momentum 
on a 12 4 (16 4 ) lattice is formed by about 0.12% (0.03%) 
of the lowest eigenvalues and eigenfunctions of the F-P 
operator. 

For the second lowest momentum qi the contribution 
of even 200 eigenmodes is far from being sufficient to 
approximate the propagator. 



V. THE PROBLEM OF EXCEPTIONAL 
CONFIGURATIONS 

We turn now to a peculiarity of the ghost propagator 
at larger (3 of which we have reported in . It was also 
seen by two of us in an earlier SU(2) study 0- While 
inspecting our data we found, though rarely, that there 
are exceptionally large values in the Monte Carlo (MC) 
time histories of the ghost propagator at lowest momen- 
tum. Those values are not equally distributed around 
the average value, but rather are significantly larger. For 
details we refer to reference 0. 

We have tried to find a correlation of such exception- 
ally large values in the history of the ghost propagator 
with other quantities measured in our simulations. For 
example we have checked whether there is a direct corre- 
lation between the values of the ghost propagator G(k) 
as they appear in the MC time histories (see e.g. Fig. 5 
in p|) and the lowest eigenvalue Ai of the F-P operator. 



In Fig. [7] we show such correlation in a scatter plot for 
different lattice sizes at p = 5.8 and 6.2. There each en- 
try corresponds to a pair [Ai, G(k)] measured on a given 
gauge copy of our sets of f c and be copies. It is visible in 
that figure, gauge copies giving rise to an extremely large 
MC value for the ghost propagator are those with very 
low values for Ai. This holds for the 12 4 and 16 4 lattice. 
However, a very low eigenvalue is not sufficient to ob- 
tain large MC values for G{k) as can be seen in the same 
figure. It is not excluded that such gauge copies with 
extremely small eigenvalues would turn out to be excep- 
tional for another realization of lowest momentum q(k) 
than those two we have used. This might explain why 
some configurations with extremely small lowest eigenval- 
ues were not found to be exceptional with respect to the 
ghost propagator at k = (1, 0, 0, 0) and k — (0, 1, 0, 0). 

In the light of Eq. it is not adequate to concentrate 
just on the lowest eigenvalues. Instead, one can monitor 
the contribution of a certain number of eigenvalues A^ and 
eigenmodes $i(/c) to the ghost propagator at some mo- 
mentum in question. Therefore, we have compared the 
truncated sums G(k\n) according to Eq. Q with the MC 
history values of the full ghost propagator G. In fact, we 
show in the scatter plots in Fig.[S]the ratios G(k\n)/G(k) 
versus G(k) for n — 10 and for various lattice sizes. 

Obviously there is a strong correlation between the 
chosen group of low-lying modes and the MC time his- 
tory values of the full ghost propagator. Indeed, if we 
consider values G{k) > 15 to be exceptional in the left- 
most panel (12 4 lattice) we find that the contribution of 
the 10 lowest modes amounts to more than 75% of the ac- 
tual value of the ghost propagator. On the opposite, for 
low G{k) values the main contributions come necessarily 



7 



12 4 16 4 24 4 



O 0.6 
£ 0.4 







be • 


• 




fc + 


• * 


- 






























-f 





6 12 18 24 10 15 20 30 40 50 

G(k) G{k) G{k) 

FIG. 8: Scatter plot of the relative contribution of the trun- 
cated sums G(k\n) over eigenmodes (see Eq. to the full 
ghost propagator values G(k) versus G(k) for lattice momenta 
k = (1, 0, 0, 0) and k = (0, 1, 0, 0). From left to right the lat- 
tice sizes are 12 4 , 16 4 and 24 4 all for (5 — 6.2. Data for fc 
and be gauge copies have been plotted separately. 

from the higher eigenmodes, while the 10 lowest modes 
contribute a minor part only. A similar but less domi- 
nant contribution of the 10 lowest modes is found for the 
time histories produced on larger lattices (16 4 , 24 4 ). 

VI. CONCLUSIONS 

In this paper we have investigated the spectral prop- 
erties of the F-P operator and their relation to the ghost 
propagator in SU(3) Landau gauge. The configurations 
under examination have been generated on a 24 4 lattice 
at (3 = 5.8 and on 12 4 , 16 4 and 24 4 lattices at (3 = 6.2. 

As expected from Ref. Q we have found that the low- 
lying eigenvalues are shifted towards A — as the volume 
is increased. The result is an eigenvalue density p{\) be- 
coming steeper rising close to A = 0. We have also shown 
that the Gribov ambiguity is reflected in the low-lying 
eigenvalue spectrum. In fact, the low-lying eigenvalues 



extracted on be gauge copies are larger on average than 
those on f c copies. Thus better gauge- fixing (in terms 
of the gauge functional) inhibits the above-mentioned 
tendency, keeping the gauge-fixed configuration slightly 
away from the Gribov horizon. 

The study of the ghost propagator in terms of the 
eigenvalues and eigenmodes of the F-P operator reveals 
that there is a dominance of the low-lying part of the 
spectrum at lowest momentum. About 200 eigenmodes 
are sufficient to reconstruct the asymptotic result up to 
a few percent at the lowest momentum on a 12 4 lattice 
at (3 = 6.2. With respect to the whole set of 8(12 4 - 1) 
non-trivial eigenvalues, this is a fraction of about 0.12%. 
For larger volumes the number of necessary eigenmodes 
seems to be somewhat larger. For the next higher mo- 
mentum, saturation needs a much bigger part of the low- 
lying spectrum. 

On average the F-P eigenmodes are not localized, how- 
ever, few large values have been seen among the lowest 
eigenmodes which so far could not be correlated to other 
measured quantities. 

Analogously with observations made in Ref. Q we have 
reported recently [5j on exceptionally large values in the 
Monte Carlo history of the ghost propagator. These we 
have seen at (3 — 6.2 and only for some lattice momenta 
k realizing the lowest physical momentum q(k). In the 
study at hand we have shown that these outliers can 
be assigned to the contribution of the ten lowest F-P 
eigenmodes to the ghost propagator at this particular k. 
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